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Abstract 

We reconsider the question of the spectral statistics of the A;-body 
random-interaction model, investigated recently by Benet, Rupp, and 
Weidenmiiller, who concluded that the spectral statistics are Poisso- 
nian. The binary-correlation method that these authors used involves 
formal manipulations of divergent series. We argue that Borel sum- 
mation does not suffice to define these divergent series without further 
(arbitrary) regularization, and that this constitutes a significant gap 
in the demonstration of Poissonian statistics. Our conclusion is that 
the spectral statistics of the A:-body random-interaction model remains 
an open question. 

PACS numbers: 02.50.Ey, 05.45. -Kb, 21.10.-k, 24.60.Lz, 72.80.Ng 
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1 Introduction 



In recent work, Benet, Rupp, and Weidenmiiller (hereafter BRW) fT], § have 
considered the spectral statistics of fc-body random interaction models P, 
^ 1^. These are models of fermions (usually; boson versions can also be 
studied 0) in which there are £ single-particle states occupied by a total of 
m particles. The models are specified by hamiltonians of the form 

Hk = ^ ^jl-jk,h:.ik^]i ■ ■ ■ ^Ik'^ik ■ ■ - (^ii ) (1) 

i<ii<i2<- •<Jfc<<? 

I<il<i2<...<ii^<e 

where the a^'s and a's are creation and annihilation operators, and the Vs 
are (for the case of the unitary ensemble) complex numbers that obey 

Vjl...jk,il--ik = ^ii...ik,ji...jk ■ 

Each independent component is a Gaussian random variable with mean zero 
and variance v^f,, 

Vh...jkM-ik^j[-3',,,i[-i'^ = ^Ik ■ ■ ■ ^jki'Jhji ■ ■ ■ Kj'k ! i^) 

the overbar denotes averaging over the ensemble, and VQk is a normalization 
constant that sets the energy scale. In this paper, for pedagogical simplicity, 
we will restrict our attention to the unitary ensemble. 

These models can be viewed as caricatures of complex systems of inter- 
acting particles, such as nuclei, multi-electron atoms, or quantum dots. The 
most physically interesting case is = 2: H2 can be thought of as represent- 
ing two-body interactions among a set of m particles occupying i degenerate 
single-particle states. 

In their work, BRW analyzed the spectral statistics of these models, and 
found that they are Poissonian for i ^ m ^ k. This is something of a 
surprise; numerical simulations of the k = 2 model have generally found 
Wigner-Dyson spectral statistics (see, e.g., [0, Even more surprising, the 
proof of this offered by BRW extends trivially to the case of nondegenerate 
single-particle levels specified hj H = Hi + H2, for arbitrary values of f 02/^^01 
i- 

This is counter-intuitive on at least two levels. First, mathematically, one 
would naturally expect that determining the spectral statistics of Hi + H2 
would be a much harder problem than it is for either Hi or H2 alone. Second, 
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physically, we expect this model to capture the essential physics of a chaotic 
many-body system (with two-body interactions), just as a single random 
matrix is often a good model of a chaotic few-body system. According to 
the Bohigas-Giannoni-Schmit conjecture |]10| (see also a classically 



chaotic system should exhibit Wigner-Dyson statistics at sufficiently high 
energy (irrespective of the number of degrees of freedom). If the BRW result 
is correct, we must then conclude that either the BGS conjecture is false for 
many-body systems, or that the two-body random-interaction model does 
not correctly capture this aspect of the physics of these systems. 

In this paper we will show that there is in fact a gap in the BRW proof, 
and that this gap is not easily bridged. Thus, we will argue, Wigner-Dyson 
statistics for these models is still an open possibility. 

In section 2 we review and expand the analysis of BRW [|l], , and explain 
how their proof of Poissonian spectral statistics might fail to hold. In sec- 
tion 3 we specialize to the case k = 1 (where the detailed analysis simplifies) 
and show that, within the context of the binary correlation method used 
by BRW, the proof does break down in the manner suggested in section 2. 
However, the binary correlation method is not exact, and the appropriate 
conclusion is that one must go beyond this approximation to obtain a reli- 
able evaluation of the spectral statistics of these models. We elaborate on 
this further in Section 4. 



2 The binary correlation method 

We begin with the resolvent 

G{z) = tr , (4) 

where tr denotes the normalized trace (tr 1 = 1) over states of m particles. 
BRW compute the connected correlation function 



Rizi, Z2) = G{z,)G{z2) - G{z,) G{z2) . (5) 
Then, since the density of states p{E) is given by 

p{E) = — \G{E - ie) - G{E + ie)] , (6) 
2iTi L -I 
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where E is real and e is a positive infinitesimal, we have 



(2vr) 



-R{zt,z^) - R{z^ ,z^ 



(7) 



where zf = Ei iie. Thus, we can extract the connected density-density cor- 
relation (which provides information on the spectral statistics) from -R(zi, ^2), 
provided we can evaluate R{zi, Z2) with both zi and Z2 on either side of the 
real axis. 

To compute R{zi, Z2), BRW expand in powers of H, 
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G{z^)Giz2) = E E ^TT^ tr(i70 tr(//' 



(8) 



r=0 s=0 ^1 



The ensemble average vanishes unless r + s is even, and, in this case. Wick's 
theorem can be used to express it as a sum over {r + s — l)\ \ pairwise contrac- 
tions of V^s. Each contraction that occurs inside one of the two traces can be 
shown [0 to yield a factor of VQj^Ak, where is a calculable number (given 
below), provided the total number of these contractions is much less than 
i. This is essence of the binary correlation method, introduced by Mon and 
French (see also [jl2|)- It is convenient to choose Vok so that fgfc^fc ~ 1' 
then we have 



A. 



m ' 



m + k 
k 



(9) 



It remains to count the total number of contractions within each trace, and 
to evaluate the contractions across the two traces. Performing the first task 
yields i 



where 



and 



R{Zl,Z2) = J29n{Zl)gniz2)Tn 
n=l 

(2p- 1)!! (n + 2p\ 



9n[Z) 



E 

p=0 



z 



2p+n+l 



n 



(10) 

(11) 



T„ = tr(i7-)tr(iJ") ; (12) 

the double overbar means that all contractions of V^s are to involve one V 
from each of the two traces. Eqs. ([T0HT2D are equivalent to eq. (68) of 0. 
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BRW show that T„ ~ ioi I ^ m ^ k, and so conclude that -R(-2i, Z2) 
vanishes in the £ — ^ 00 hmit; this is indicative of Poissonian spectral statis- 
tics. However, this conclusion is suspect if we do not first understand the 
convergence properties of the series in eq. ([T0|). To see why this is necessary, 
consider the mathematical example 



^ n\ 



^=E^, (13) 



n=l 



where £ is a positive integer. Clearly, in the limit £ — > cxd, each term of this 
series vanishes, and so it is tempting to conclude that S' = 0. However, we 
can evaluate S by Borel summation. We use 



n! = / dte-'e (14) 
in eq. (|l^) and do the sum to get 



S= / dte-'- . (15) 

7o ^-t ^ ' 

This integral is not defined for positive real I. We could attempt to define it 
by analytic continuation, but there would still be an ambiguity, corresponding 
to whether the positive real axis is approached from above or below. More 
importantly, in our case £ is a positive integer that counts the number of 
single particle states; therefore it does not seem to make sense to analytically 
continue to complex Our conclusion in this case would be that S is simply 
not defined by the Borel procedure. We wish to examine whether the same 
problem arises for the series in eq. (jT^). 

The series for gn{z) in eq. ([TT|), on the other hand, can be defined for any 
z with Im^ 7^ by Borel summation followed by analytic continuation. To 
demonstrate this, we begin with the combinatoric identities 

'" + 2''>' and (2p-l)..= <M. (16) 




and so 



n!(2p)! ' ' 2Pp\ 



^One could raise the same objection to dimensional regularization in quantum field 
theory, but there the results can be verified by a variety of different and more physically 
motivated schemes. 
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We now use eq. (|T^) with n n + 2p to get 



p: 



c/te-*re*'/2-' . (18) 



The integral converges provided R,ez~^ < 0. This condition is satisfied for 
z = ±i\z\e^'^ with — ^vr < < +^7i. We now rotate the t contour in the 
complex t plane so that it runs along a straight line from zero to complex 
infinity at an angle of relative to the positive real axis. Then we set t = re^'^, 
where r is real and runs from zero to infinity. We now have 



where the phase in the numerator of the prefactor comes from the change of 
variable t —>■ re**^, and the phase in the denominator comes from z = ±i\z\e^'^. 
We now change the integration variable to u = t/\z\ to get 

\n+l ^00 , . , „ , 



9n{z) = 



n 



n+l 



/ rfMe^^""M"e-" (20) 
Jo 



This integral converges for all z, and so constitutes an analytic continuation 
of eq. (|l8|); the ± symbol should be interpreted as the sign of Imz. Thus 
gn{,z) is discontinuous across the real axis, and if we take z = z^ = E ±ie, 
we get 

g^{z^) = ^^jp due^^'^-u-e--'/' . (21) 
We also note that, from our eq. (|ll|) and eq. (65) of 0, 



G{z)=g,{z). (22) 
Then, from eq. (|^) and eq. (0), we have 
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p{E) = — 9o{z )-9o{z^ 



fOO „ 

— / du (e-*^" + e+'^") e-" /2 
2tt Jo 

^ :exp(-EV2) . (23) 



This is the classic result of Mon and French : the ensemble-averaged den- 
sity of states of the /c-body random-interaction model is Gaussian. 
Returning to R{zi, Z2), we use eqs. (|10|) and (^) to get 



R(z+ 4) = ^ / dudv e^(^i-±^2.) e-(" +^ , (24) 

Jo 

where we have defined 

00 n 

ny)^Y.^j^.T^- (25) 

We see that, in order for both R{zi,Z2) and R{zi,Z2) to be well defined, 
F{y) must be free of singularities on the real axis (positive or negative). Only 
then can we perform the integrals over u and v in eq. (p^ for both R{zf , Z2) 
and R{zf,Z2) without further (arbitrary) regularization. 

To see whether or not this obstacle arises, we must evaluate T„. Here we 
have an immediate difficulty. We have already invoked the binary correlation 
approximation in eq. (pJ^; the terms in this series receive corrections when the 
summation index p becomes comparable to i, and there is no straightforward 
way to calculate these corrections exactly. Similarly, in section 3 we will 
evaluate T„ (for k = 1), but our method will require n <^ i. It is therefore 
unsuitable for reliably determining the asymptotic behavior of the series in 
eq. (pSl). Still, it is worthwhile to see whether or not the problem of a 
singularity on the real axis arises within this approximation. We therefore 
turn to the calculation of T„ for k = 1. 

3 Analysis for k = 1 

We specialize to the case k = 1 (and drop the corresponding "1" subscripts): 

H = J2Vjia]ai, (26) 

with 

VjiVj'i' = vl SjiiSiji . (27) 

Of course, we already know the answer for this case: the spectral statistics 
are Poissonian, because the spectrum simply consists of the linear sum of 
the m single-particle energies that are obtained by diagonalizing the ^ x 
£ hermitian matrix V . These single-particle energies obey Wigner-Dyson 
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statistics, but their sums (for m ^ 1) obey Poisson statistics. However, 
the binary correlation method can still be used, and it is important to see 
whether or not it gives the correct answer (or any answer at all). 

We first introduce some shorthand notation. Let / = and 

J = {ji,---,in}; let 

5i,j = . . . . (28) 

Let P, Q, and R denote permutations of the n indices of /. Then, from 
eq. (H), we can write 

(^.n---^„.J(%;...l^;.J = vl- Si,Rj'hRi' (29) 

R 

Another useful bit of shorthand is 



Thus eq. (|l^) becomes 

Tn = t'o'" E E Si,RJ'h,Rr triA^Aj) tTiA\,Aj,) . (31) 

R IJI'J' 

The trace is a f/ {£) invariant operation, and so we must have 

ti{AlAj)=J2CpSj,pi (32) 
p 

for some set of coefficients Cp. Substituting this expansion into eq. (0), we 
get 

Tn = E CpCq E ^I,RJ'^J,Rr^J,PI^J',QI' 

P,Q,R IJI'J' 

= vl- E CpC'QE^/.iiQRP/ (33) 

where R is the inverse of R. 

The expression X)/ h,Pi can be evaluated in terms of the cycles of P. 
Consider, for example, the permutation 2431 of 1234; it contains one cycle of 
length one (since 3 remains in its original position), and one cycle of length 
three (since 1 is replaced by 2, 2 is replaced by 4, and 4 is replaced by 1). 
Let 7c(-P) denote the number of cycles of length c in permutation P; for 
our example, 71 = 73 = 1 and 72 = 74 = 0. Since each element of any 



8 



permutation is in exactly one cycle, we have J2c=i1c{P)c = n. In J2i^i,pi, 
each cycle of P ultimately results in a factor of J2i = ^- Thus we have 

Y^Sj^Pj = e^MP) , (34) 

and hence _ 

Tn = vr E CpCQi^M^^^P) . (35) 

P,Q,R 

So far we have made no approximations in our evaluation of T„. We now 
notice that, for £ ^ 1, the dominant term on the right-hand side of eq. (pS]) 
is the one with the largest value of J^c'lciRQRP)- This occurs when 71 = n 
and 7c = for 2 < c < n, which in turn implies RQRP = I, or Q = RPR. 
Thus we have 

T^^vl^rY^CpCj^Pn. (36) 
pfi 

where = denotes equality in the limit of large £. 

Next we must evaluate Cp. Starting with eq. 
sum over I to get 

Etr(4,A,) = E^pEW 
/ PI 

p 

= Cqt . (37) 

To get the last line, we used the same large-£ argument that gave us eq. (^6|). 

Now we need to evaluate X)/ ti'(^p/^/)- To do so, imagine rearranging 
the a^'s so that their indices are in the same order as the a's, and interleaved 
among them in the standard pattern of eq. (PD|). Since aja] = — a]a|, we see 
that, before accounting for the presence of the a's, rearranging the a^'s in this 
way will yield a factor of (—1)^. Passing a^'s through a's yields Kronecker 
deltas, since a^a] = — (a]aj — (5jj), but no possibility of an extra overall minus 
sign, since the total number of exchanges of a^'s with a's is always even. Let 
us temporarily ignore the Kronecker deltas. Then, once all the a^'s have been 
put next to their partner a's in the standard pattern, we have X^i A'^i = ^■ 
Thus, ignoring the Kronecker deltas generated by moving a''"'s past a's, we 
have Cp = (—l)^m"'. 

For n <^ i, most of the Kronecker delta's can indeed be neglected; the 
exception occurs when one of them has two identical indices, leading to 



32), we set J = QI and 
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J2i ^ii = ^- To account for these dominant Kronecker deltas, we first note 
that, for a given permutation P, 7i(-P) is the number of elements that are 
left in the same place by P, and these contribute no Kronecker deltas. Let 
r{P) be the number of elements that are moved to the right of their original 
positions by P; only these a^'s need to pass through their partner a's on 
being returned to the standard order. For these, the factor of m must be 
replaced hj m — i. Finally, the number of elements that are moved to the 
left of their original positions by P is n — 71 (P) — t(P), and these contribute 
no large Kronecker deltas on being returned to the standard order. Thus we 
have 

Cp ^ {-if m^l(^) ^n-7i(P)-r{P) _ ^yiP) ^gg^, 

for n <^ i. 

Next we notice that the inverse permutation P has the counts of the left- 
and right-moving elements exchanged. Thus, 

Cp ^ (-1)^ r" m^i(^) m'-(^) (m - , (39) 

Then, using the fact that (—1)^ and 71 (P) are each invariant under unitary 
transformations of P, we have 

Putting eqs. (53) and (BO) into eq. (PB|), we get 



m \Jl(P)-r{P)+r{RPR) 

Tn = vr[m{m-i)rr-Y. (-—-.) ■ (41) 

The convention of eq. (P) with k = \ yields v^\m{m — £)]" = (—1)". Then 
we can write 

T„/(n!)2 = r"M„(//(l-/)) , (42) 
where f = m/l is the filling fraction, and we have defined the function 

M^{X) ^ ^ Y^^_^y.iP)-riP)+r(nPR) (43) 

{n.) pp 

Remarkably, it is possible to evaluate M„(a;) exactly for arbitrary n without 
performing the sum over the permutations P and R in eq. (^31) explicitly. 
This calculation is given in the Appendix. We find that Mn{x) is everywhere 
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positive and convex; it diverges as x — and x +00, and has a minimum 
between x = and x = 1. While we are not able to evaluate the asymptotic 
form of Mn{x) analytically, numerical evaluation reveals that, for large n, 

Mnix) ~ nA" , (44) 

where A is a number that depends on x; if x is real and positive, so is A. 
Using eqs. (|4^) and (H ) in eq. (p5|), we find 



Fiv) ^^/^ (45) 

which obviously has a singularity on the positive real axis at y = +£/X. 
Hence the integral in eq. (^) for R{zf,Z2) is not well defined. Therefore, 
we argue, the binary correlation method fails to give us the density-density 
correlation function. 



4 Conclusions 

It is significant that, while we cannot compute R{zi , Z2) by the binary cor- 
relation method, we can compute R{zf , ^2")- In this case (at least for k = 1), 
there is no singularity in the integrand in eq. (P^ ) , and so we can complete the 
integrals, take the i ^ 00 limit, and conclude that R^zf, Z2) = R{zi , Z2Y = 
0. This is consistent with other approaches to the computation of this spec- 
tral correlation function. For example, in the periodic-orbit approach to 
chaotic systems, one finds (see, e.g., [[l^] ) 



G{z^) =Y.Wp e^'^^/^'e-'^- , (46) 

p.o 



where Sp is the action of the orbit (including the contribution of the Maslov 
phase), Tp is its period, and Wp is the Gutzwiller weight factor. The key 
point is that, in the limit of small G{zi)G{z2) contains only large (^ 27r) 
phases, and so summing the orbits should give zero. On the other hand, 
G{zi)G{z2) has negative phases from the first factor and positive from the 
second, so cancellations can occur; in the most straightforward approach (the 
so-called diagonal approximation), one simply keeps only those terms with an 
exact cancellation. This is adequate for computing the density-density corre- 
lation function at large l-Ei— for chaotic systems, it is nonzero, and agrees 
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with the Wigner-Dyson prediction of random-matrix theory. Thus, for the 
k-hodj random-interaction model, the binary-correlation method correctly 
and unambiguously gives us R{zi,Z2) = 0, but, as we have argued, formal 
manipulations of divergent series are not adequate for extending this result 
to R{zf,Z2)- We have seen explicitly how the binary-correlation approxi- 
mation fails to give us a series for R{zf, Z2 ) that can be evaluated by Borel 
summation in the case k = 1. 

It is important to remember, though, that the binary correlation approx- 
imation breaks down precisely where we wish to apply it: when the number 
of contractions approaches the number of single-particle levels. Thus, to 
reach a firm conclusion one way or the other, we are forced to go beyond this 
approximation. At the present time, doing so appears to present a severe 
challenge to the available technologies. 

One of these technologies, also investigated by BRW, is the supersymmet- 
ric sigma model BRW found that, at the tree and one-loop levels, the 



sigma model predicts Wigner-Dyson statistics for A; <^ m ^ £, with correc- 
tions that go to zero as Unfortunately, this does not settle the matter, 
because this prediction also applies to = 1, where it is known to be false. 
One must therefore presume that higher- loop corrections are important. 

Thus the only conclusion that seems completely safe at this juncture 
is a disappointing one: the nature of the spectral statistics of the /c-body 
random-interaction model remains an unsolved problem. 
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5 Appendix 

We wish to evaluate 

{-lY{n\fM^{-x) = J2 ^ (47) 

P,R 
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where P and R are permutations of n elements, 71 (-P) is the number of 
elements left in their original positions by P, and r{P) is the number of 
elements that are moved to the right of their original positions by P. For 
small values of n, we can perform this sum over (n!)^ terms explicitly. We are, 
however, interested in its behavior for large n, and so this strategy quickly 
becomes untenable. 

Let us warm up by computing 

K(a;)=5:x^'(^) . (48) 
p 

Let Nn^r be the number of permutations of n elements in which exactly r of 
the elements move to the right. Then 

n-l 

Rnix) = ^iV^X • (49) 

r=0 

We can compute Nn^r by the method of inclusion and exclusion Let 
Sn,r be the number of permutations in which at least r elements move to the 
right. Then 




Nn,r = T.i-^)'\ T ]Sn,r+k. (50) 



To obtain Sn,r, let us first count the number of permutations in which r 
specific elements move to the right, irrespective of what happens to the re- 
maining elements. Actually, since we conventionally count from left to right, 
it is slightly more calculationally convenient to count the number of permu- 
tations that move r specific elements to the left, rather than to the right; 
clearly, by left-right symmetry, the choice of direction is irrelevant to the 
count. So, let A„ be the number of permutations in which the if^ element 
moves to the left. There are ii — 1 possible places for the if^ element to go, 
and (n — 1)! permutations of the remaining elements; thus 

^n,n = («i-l)(^-l)! . (51) 

Similarly, let An,ii...iri with ii < . . . < be the number of permutations in 
which all the named elements move to the left; by a similar argument, we 
have 

An,r,..A. = {h - 1)(^2 - 2) . . . (v - r)(n - r)! . (52) 
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Now Sn,r is given by the sum over all possible values of the i's of An^i^,„i/- 

n 43 — 112 — 1 

Sn,r = 5^ • • • An,ii...ir 
ir=r 42=2 ii=l 

n 13 — 1 22 — 1 

= (n-r)!E(v-r)...5:(z,-2)E(zi-l) 

ir=r «2=2 *i=l 

n—r js j2 

= {n-ry.Y. Jr---J2 ^ (53) 

jr=0 J2=0 jl=0 

where, in the last line, ja = ia — a- Eqs. (^), (0), and ( ^31) give us Rn{x). 
The first few of these polynomials are 

Ri{x) = 1 
Riix) = 1+x 
R^(x) = l + 4x + x^ 

R^(x) = 1 + 26x + 66x2 + 26x^ + . ^^^^^ 

Also, it turns out to be convenient later if we adopt the convention 

Ro{x) = l/x . (55) 

Now consider J2r x^'^^^^^ for fixed P. Let P belong to cycle class C, 
specified by the numbers 7c of cycles of length c, with Ylc=i ^Ic = The 
permutation RPR belongs to the same class (since class is preserved by uni- 
tary transformations), and furthermore each permutation in the class appears 
equally often as R is varied over the group. Thus 

,r(RPR) _ ^r{P) 



^^.(HP^)^_1,^^.(P), (56) 

where 



R p^c 



'^^ ^ nOT 



is the number of permutations in class C |T5[. Let Nc^r be the number of 
permutations of the elements in a particular cycle of length c in which exactly 
r of the elements in this cycle move to the right. A cycle can be specified 
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by an ordering of its elements with the smallest first; thus a cycle of length 
c can be mapped to a permutation of the remaining c — 1 elements. This 
mapping can be shown to imply Nc^r = ^c-i,r-i- Thus, 



c-l _ 
r=l 



.X 



c-l 

Yl ^c-i,r-ia;' 

r=l 

c-2 

s=0 

xRc-lix 



s+1 



(58) 



The convention of eq. (^) correctly treats the case c = 1. Also, the number 
of ways of assigning elements to the cycles is 



Na 



Thus 



P&C 



r{P) 



n"=i(c!)>7c! ■ 
Na n xRc-i{x) 

cycles 
n 

NAl[[^Rc-iix)Y' 
xRc-lix 



(59) 



c=l 

„ ^ 



c=l /c- 



C! 



7c 



(60) 



Therefore, by eqs. (|56[) , (|57[) , and (^ 



E 



r(RPR) 



c=l 



x-Rc_i(x) 



7= 



(61) 



Now we can write 

(-l)"(n!)2M„(-x) 



P,R 

c Pec 



,71 (P)-r(P)+r(iJPR) 



71 (P)-r(P) 



E 

R 



X 



r{RPR) 



(62) 



The sum over R is given by eq. (|6TD , and the sum over P G C is given by 
eq. ( |60|) with x-Rc-i(a;) replaced by x'''=i~^i?c--i(l/a;)- The sum over classes 
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is equivalent to a sum over all possible values of each 7c with the constraint 
that J2c=i C7c = n. This is most easily treated via a generating function 



n=0 



Mix,z) = Y.Mn{x)z 

Then, from eqs. (|60H 6^) we have 
M{x, z) 



(63) 



OO OO -| 

7c=0c=l 'c- 
.c=l 



'-zY{-xY^^R,_A-x)R,. 



-l/x) 



exp 



c!(c- 1)! 
-xY^^R^_r{-x)R,^i{-l/x) 



die- I] 



(64) 



The expansion of eq. (6|) in powers of z yields M„(x) as the coefficient of 
z^. The results for small n are 



Mi(x) 
M2(x) 
M^{x) 
Mi{x) 



X 



(l + x2)/2 

(l/a;-2 + 7x + 2x^)/(2!3!) 

(1/x^ - 8/x + 48 - 32x + 49x^ + 6x^)/(3!4!) . 



(65) 



We have verified this procedure for computing M„(x) up through n = 7 by 
comparison with the brute-force summation of the right-hand side of eq. (^T]). 

It is tempting to conclude from eq. ( |65D that, in the limit of small x, we 
can take M„(x) = a;~*-"~^V(^K^~l)0- Small x corresponds to small filling 
fraction / = m/£, which is what we want if we are interested in the limit 
£ — i> OO with m fixed. This limiting form for M„(x) would lead to convergence 
of the series in eq. (|25|) for all ?/, with no singularities on the real axis, thus 
apparently validating the BRW result. However, including the first correction 
yields 



MJx) 



X 



-(n-2) 



/{n\{n-l)\ 



(66) 



1 - iT-2n)x + 

and so we see that we must have x <^ 2^" before the first term alone is 
an adequate approximation. Instead, for the purpose of determining the 
convergence of the series in eq. (P5|), we should take the limit of large n with 
X held fixed. We have not found a way to do this analytically, and so have 
resorted to numerical methods. For a given value of x, we can evaluate Mn{x) 
up to around n = 500 from eq. (0) in a reasonable amount of computation 
time. We find, for x = 0.001, 0.01, 0.1, 0.5, and 1.0, that the behavior of 
Mnix) for > 10 is very well fit by eq. (||), with A = 17.4, 3.28, 0.797, 
0.435, and 0.405, respectively. 



16 



References 

[1] L. Benet, T. Rupp, H.A. Weidenmiiller, |cond-mat/0010425| , Phys. Rev. 
Lett. 87, 010601 (2001). 

L. Benet, T. Rupp, H.A. Weidenmiiller, |cond-mat / 00 1042B| , Ann. Phys. 
(N.Y.) 292, 67 (2001). 

J.B. French and S.S.M. Wong, Phys. Lett. B 35, 5 (1971). 

O. Bohigas and J. Flores, Phys. Lett. B 34, 261 (1971). 

K.K. Mon and J.B. French, Ann. Phys. (N.Y.) 95, 90 (1975). 

T. Asaga, L. Benet, T. Rupp, H.A. Weidenmiiller, |cond-mat /010736"3| , 
Europhys. Lett. 56, 340 (2001); |cond-mat/0107364| , Ann. Phys. (N.Y.) 
298, 229 (2002). 

T.A. Brody, J. Flores, J.B. French, P.A. Mello, A. Pandey, and S.S.M. 
Wong, Rev. Mod. Phys. 53, 385 (1981). 



T. Guhr, A. Miiller-Groeling, and H.A. Weidenmiiller, |cond- 
mat/970730ll , Phys. Rep. 299, 189 (1998). 

L. Benet, T. Rupp, H.A. Weidenmiiller, private communication. 

O. Bohigas, M.-J. Giannoni, and C. Schmit, Phys. Rev. Lett. 52, 1 
(1984). 

G. Casati, F. Valz-Gris, and L Guarnieri, Lett. N. Cim. 28, 279 (1980). 

J.J.M. Verbaarschot and M.R. Zirnbauer, Ann. Phys. (N.Y.) 158, 78 
(1984). 

E. Bogomolny, |nlin. CD/00 11023| , in New directions in quantum chaos: 
proceedings of the International School of Physics "Enrico Fermi", 
course CXLIII, G. Casati, L Guarneri, and U. Smilansky, eds. (lOS 
Press, Amsterdam, 2000). 

K.B. Efetov, Adv. Phys. 32, 53 (1983). 

J. Riordan, An introduction to combinatorial analysis (Wiley, New York, 
1958). 



17 



